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A method of determination of the Gibbs free surface energy and the radial 
dependence of adhesion forces of adsorbents in respect to the interfacial 
water was developed on the basis of measurements of the 'H NMR signal 
intensity of unfrozen water at T « 273K. It was considered features of 
application of this method to highly disperse oxides, modified oxides, 
microporous and mesoporous adsorbents. The free surface energy values 
are computed for a variety of oxide and carbon-mineral materials. 


1. INTRODUCTION 


NMR spectroscopy is widely used to study the structure of solids and 
asignificant portion of the results has been obtained utilizing the solid-state 
spectroscopy of high resolution with magic-angle spinning and 
cross-polarization (CP-MAS) or combination of the rotation with the 
multi-pulse spectroscopy (CRAMS). Exploration of interaction between 
adsorbed molecules and solid surfaces can be based on the measurements of the 
temperature dependencies of longitudinal and cross-relaxation of nuclear spins 
of adsorbates [1]. The mobility of adsorbed molecules decreases with 
strengthening adsorption interactions, which reduce the relaxation time. 
Intermolecular (adsorption) interactions can be described in details on study of 
the self-diffusion measured using pulse-gradient of the magnetic field [2-5]. 
These measurements demonstrate that interaction between adsorbate and 
adsorbent is propagated through a relatively thick layer, which can reach 
2-3 nm or larger [3,6]. In aqueous suspensions of oxides, interaction of water 
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with solid surfaces can be detected even at 100 nm distance due to long-range 
component of intermolecular interaction [7], formation of the electrical double 
layer [8-10], and electrostatic and polarization effects [11,12]. 

An important parameter characterizing interaction between solids and 
liquids is the Gibbs free surface energy (ys); however, its exact measurement is 
difficult [13-20]. In liquid water 


Ys = [AG(x)dx (1) 


(where x is the thickness of the bound water layer, AG is the change in the 
Gibbs free energy dependent on x) equals to summarized changes in the Gibbs 
free energy of solids/water due to the existence of the interfaces. The 'H NMR 
spectra of adsorbed water are useful and informative on the study of interfacial 
water, which can be frozen at T < 273K due to the impact of solid surfaces 
reducing the free energy of the water. The layer thickness of water unfrozen at 
T < 273K can be estimated using the 'H NMR signal intensity measured during 
freezing and de-freezing processes. The interfacial water is in a quasi-liquid 
state at T < 273 K and characterized by a relatively narrow 'H NMR signal. The 
signal of ice (as well as OH groups of solids) is not observed in the 'H NMR 
spectra because a large difference in the times of cross-relaxation of protons in 
liquid and frozen waters. 

Dependence of the unfrozen water layer thickness (i.e. the intensity / of the 
'H NMR signal) on temperature in aqueous suspensions or hydrated powders 
was studied by many authors [21-32]. It was shown that this thickness depends 
on the hydrophilic properties of adsorbent surfaces; for instance, it decreases as 
follows silica gel > kaolinite > Teflon [28]. There is a /(T) hysteresis loop on 
freezing-de-freezing [23,31], which can appear due to overcooling of water on 
freezing [23]. Therefore, the /(T) functions obtained on de-freezing of samples 
after freezing to T lower than a minimal temperature of water overcooling are 
more appropriate to determine the characteristics of the interfacial water. 
Computation of changes in the Gibbs free energy using the dependence of the 
signal of unfrozen water on temperature is based on the fact that the freezing 
temperature of the interfacial water depends on the free surface energy of 
adsorbent and the distance between water molecules and solid surfaces [32-35]. 
Condition of water freezing at the interface is the equality of the free energies 
of water and ice, and typically the closer the adsorbed molecules to the 
surfaces, the lower the temperature of water freezing. 
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2. CHARACTERISTICS OF THE BOUND WATER LAYERS 
DETERMINED USING 'H NMR 


It is of interest to consider the application of this method to the aqueous 
suspension of fumed silica (specific surface area Sper ~ 300 m’/g), previously 
pressed (~1.14x10° Torr) and ball-milled. Figure 1a shows the 'H NMR signal 
of the unfrozen water in such a suspension at different temperatures. The 
observed signal is linked to mobile water molecules affected by the silica 
surfaces. This signal is individual at the chemical shift 54 ~ 4 ppm (in respect to 
TMS), whose intensity decreases but half-width increases with lowering 
temperature in consequence of a reduction of the molecular mobility. The 
concentration of unfrozen water (C,,) can be estimated using the signal 
intensity comparing with that of a given amount of water adsorbed on this 
oxide from the gas phase. To determine the concentration of water adsorbed on 
a sample (C,), one can utilize the graph of the intensity of 'H NMR signal 
versus the amount of water (C,) added to the powder using a metering 
micropump (Figure 1b). The point of intersection of /(C,+C,) with the X axis 
corresponds to the initial concentration of adsorbed water. The thickness (d) of 
unfrozen water (as the number of the statistical monolayers) can be computed 
from the concentration of adsorbed water (C,,,) and the specific surface area of 
adsorbent 


d-3.1C,/S (2) 


Figure 2a shows the graph of the intensity of the unfrozen water signal versus 
temperature. Changes in the Gibbs free energy AG linked to lowering 
temperature of water freezing can be easily computed, as the thermodynamic 
functions of ice are known over a wide temperature range [36]. The relationship 
between AG of ice and T can be written as follows 


AG = 0.036(273 — T) (3) 


The /(T) graph for the unfrozen water can be easily transformed to the 
dependence of AG for adsorbed water on the distance of the boundary between 
frozen and unfrozen waters to the adsorbent surface (Figure 2b). 

Sometimes an increase in the intensity of the NMR signal with lowering 
temperature due to enhancement of nuclear level population (Curie rule [37]) 
should be considered. There is a relationship between the nuclear 
magnetization (M) of the system with N spins with the spin value s and the 
gyromagnetic ration g 


M = KNg*h’I(s + 1)Bo/3kT (4) 
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where B, is the magnetic field intensity, K is a constant dependent on NMR 
measurement conditions, h is Plank constant, k is Boltzman one. Notice that for 
thick water layers adsorbed on the surfaces of dispersed oxides, the choice of 
measurement conditions allows one to obtain the NMR signals of water (with 
no freezing) constant over a wide range of the temperature, i.e., the Curie rule 
can be ignored. 
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Fig. 1. (a) 'H NMR signal of unfrozen water in aqueous suspension of fumed silica, 
previously pressed and ball-milled (A,); and (b) the graph of the intensity of 'H NMR 
signal versus the amount of water (C,) added to the powder 


The capillary effects in the aqueous suspensions of fumed oxides possessing 
nonporous primary particles are absent, and AG(C,,,) determines changes in the 
Gibbs free energy of the interfacial water as a radial function of the unfrozen 
water layer thickness around spherical primary particles as C,, is proportional 
to the layer thickness x. The value of the free surface energy of adsorbents in 
the aqueous suspensions can be determined using the AG(C,,,) function, and 
Eq. (4) can be re-written as follows 


Cw 
Ys=K, JAG(C,,)dC,, (5) 
0 


where K, is a scale coefficient, C,,"^' is the thickness of the unfrozen water 
layer extrapolated at AG — O (if AG is in kJ/mol and ys in mJ/m^ that 
K, =55.6/S, where S is the specific surface area of adsorbent). Two portions of 
the AG graph can be separated (Figure 2b). The first portion corresponds to 
alarge reduction in the unfrozen water layer over a narrow AG range (near 
273 K) and the second one is linked to small changes in the thickness over 
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abroad AG range. The first type of the interfacial water can be assigned as 
weakly bound water (thickness d,) due to a long-range component of the 
surface forces. The second type corresponds to a greater change in the free 
energy (AG,) over a short distance (d,) from the surface, i.e., this is strongly 
bound water. The thickness of these layers (d, and dẹ for strongly and weakly 
bound waters) and overall changes in the free energy for them (AG, and AG,,) 
can be estimated using extrapolation of the plots to the axes. In the case of 
a nearly linear shape of these AG(d) graph portions, the free surface energy at 
the interface of oxide/water/ice can be determined using a simple equation 


Ys = Ki(AG;d, + AG,d,)/2 (6) 


This value can be computed as the sum of AG for all monolayers of the 
interfacial water 


PE Y AG,d, (7) 


where AG; is the average value for the i monolayer and 1 < i < dy = d,+d,. The 
AG, value corresponds to changes in the free energy of the first monolayer. The 
overall AG,,, value can be estimated using the AG(d) graph at 
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Fig. 2. (a) Dependence of the 'H NMR signal intensity of unfrozen water versus 
temperature; and (b) dependence of variation in the free energy of water at the 
silica/water interface as a function of distance between the surface and the 
frozen/unfrozen water boundary in the aqueous suspension of A, 


3. RELATIONSHIP BETWEEN THE FREE SURFACE ENERGY Ys 
AND OTHER THERMODYNAMIC PARAMETERS 


Immersion heat. The differential immersion heat A is the heat evolved due to 
addition of an infinite amount of liquid onto the surface at a given 
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coverage x [38]. The A value as well as AG characterizes the surface force field. 
The value A = 0 corresponds to full compensation of the surface energy by the 
energy of adsorbate interaction with the surface. The A value as well as the 
differential Gibbs free energy AG for different layers decreases farther and 
farther away from the surface. The differential Gibbs free energy equals to the 
differential work of adhesion 


AG = -W, (8) 


The Gibbs free energy increases with increasing distance from the surface and 
this increment equals to a reduction of the surface Gibbs free energy of 
immersion. The change in the Gibbs free energy of liquid is linked with the 
thermodynamic activity (a,) 


AG = -RT Ina, (9) 


According to the Gibbs-Helmholtz equation, the differential heat of immersion 
can be written as follows 


A = dAH/ox = AG — T(9AG/OT) (10) 
At the same time, the derivative of AG/T with respect to temperature is 
dAG 
acit) 7 A6*TC34) 
SECO JUPE ung. (1) 
oT I 
From Eqs. (10) and (11), one can obtain 
9 AG d 
A=-T* ——= RT? (—In 12 
oT T (gr ne. n 
The thermodynamic activity of liquid at a given coverage of the surface is 
ax = pps (13) 


where p, is the pressure of saturated vapor over the interface liquid at coverage 
x, p, is the pressure of saturated vapour over the pristine liquid. Using Eq. (10) 
one can obtain 


A = RT (ln pJ/0T), — RT (dn p/aT) (14) 
where the first term in the right part is the liquefaction heat Q, of the vapor on 


the surface of the liquid, the second one is the liquefaction heat (Qo) under 
standard conditions, then 


ÀA-2Q. - Qo (15) 
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Consequently, the differential heat of immersion equals to the difference 
between the liquefaction heats of the water vapour over the adsorbent surface 
and the bulk water. Thus, there is a similarity between the determination of the 
immersion heat in Eq. (15) and AG described above. Both the immersion heat 
and differential free energy are determined by the difference in the 
thermodynamic functions for the phase transfer to free and interfacial liquids. 
Since the free energy and enthalpy are linked by equation 


AG = À - TAS (16) 


which can be re-written for overall values of the free surface energy and 


summary heat of immersion (A;) with consideration of changes in the 
thermodynamic functions over the total thickness of the bound water layer 


Ys = A; — [TAS (x)dx (17) 
and 
À = [Adx (18) 


The dependence of the heat effect on temperature can contribute Eq. (17). The 
overall value of the immersion heat is determined at a fixed temperature; 
however, the measurement of the free surface energy by means of freezing-out 
of the bulk of water is performed over a broad temperature range. Therefore, 
Eqs. (16) and (17) can be re-written 


AG = A(xT) - TdS (19) 
Ys = JJA(x,T)dxdT — J {TdS(x,T)dxdT (20) 
Tx Tx 


In the literature, there are no data related to the dependence of the free 
surface energy of silica on temperature. However, similar investigations 
performed by means of inverse gas chromatography (IGC) for graphite [18] 
show that a dispersion component of Ys decreases from 140 mJ/m! to 125 mJ/m* 
with elevating temperature from 300 K to 360 K. One can assume that similar 
small changes in Y; can occur in the case of dispersed oxides. Then 
non-isothermicity of freezing-out of the bulk water under the 'H NMR 
measurements does not give large errors in respect to Ys, which can be 
compared with the free surface energy determined using other methods. 


Disjoining pressure and surface forces. One of frequently used 
characteristics of the surface forces is the disjoining pressure, which is the 
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difference between the pressure p, on the liquid layer located between two solid 
surfaces and the pressure pp in the bulk of liquid [7] 


Il(x) = pi — po (21) 


where x is the thickness of the liquid layer. Besides, the disjoining pressure 
equals to changes in the Gibbs free energy over the unit distance to the surface 
at constant temperature (T), pressure (p) and chemical potential () 


I(x) F= —(dG/dx)rpu (22) 


Typically the disjoining pressure is measured using direct methods to 
determine the force needed for closing two plates in liquid media and these 
methods are described in detail elsewhere [7,8]. Formally, the disjoining 
pressure value can be estimated from the data shown in Figure 2b. However, 
Eq. (22) corresponds to the determination of the II value at a constant 
temperature. Therefore for isobaric freezing-out of the bulk liquid, one can use 
an effective value of surface forces, which can be calculated using the 
relationship 


F = (AG m)/x (23) 
where m is the mass of unfrozen water in moles, x is the distance to the surface. 


If the AG value is in kJ/mol, C,,, in mg/g, x in nm, then 


AGC,,,-10° 
== —(Nig (24) 


F 


The derivative value of AG with respect to the thickness of the bound water 
layer (x) or the concentration of unfrozen water 


0 = d(AGydx ; or 0 = d(AG)/d(C,,) (25) 


Maxima of this dependence determine the characteristic distances to the surface 
(or concentration of bound water), for which changes in the type of surface 
forces are observed. For fumed silica pressed at p = 150 bar, the dependences 
0 = f(x) and F = f(x) are shown in Figure 3, where F is in GN per gram of 
adsorbent and shows the interaction force between liquid water and gram of 
silica. This value divided by the specific surface area corresponds to the 
adhesion force per unit of the surface area. In the case of direct measurements, 
the force between two plates in liquid water is estimated [39-43], but in our 
measurements, the adhesion force value is computed as the interaction force 
between the solid surface and liquid medium. These forces are the main 
characteristic of many parameters studied for composites and filled matters. 
One can separate long-range and short-range components in the adhesion forces 
for the aqueous suspension of fumed silica (Figure 3). This suggests the 
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possibility of the determination of strongly and weakly bound waters. Change 
in the type of the surface forces occurs at mean distance to the surface of 
approximately 1 nm. 


F*10E9 (N/g) 
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Fig. 3. Radial dependencies of (a) adhesion forces and (b) O(x) value for Ap aqueous 
suspension 


4. HIGHLY DISPERSE OXIDES 


4.1. Influence of the adsorbent pretreatment on adhesion forces. The 
samples of fumed silica A-300 (Pilot Plant of the Institute of Surface 
Chemistry, Kalush, Ukraine; specific surface area Sger = 300 m2/g, apparent 
density C, = 40 mg/cm’) were studied such as non-treated (A,,), heated at 520 K 
for 3h (A). wetted with water (A, ) or hexane (supplied by Aldrich) (Apex) then 
dried at 293 K (C, = 90 mg/cm’), pressed at 15 MPa, (C, = 250 mg/cm’), (A, ) 
[44,45]. 

Changes in the 1H NMR spectra as functions of the temperature and the 
silica concentration for the frozen suspension of A, are given in Figures 4a and 
4b, respectively. The spectrum of water adsorbed on A, (Figure 1a) has a single 
signal, whose width increases while the temperature goes down due to a decline 
in the molecular mobility of adsorbed water and the intensity falls as the 
thickness of region aqueous close to the silica surface decreases through 
freezing. The same type of temperature dependent changes is observed in the 
spectra registered for A, and Aex The signal of adsorbed water for A, and A, 
is observed againts the background of a considerably wider signal having 
a shape different from the Gaussian one. Because of their similarity only the 
spectra for water on A, are given. The intensity of a broad component quickly 
falls as a function of decrease in temperature, and it is not observed in the 
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spectra at T « 255 K (Figure 4). It has been shown earlier [46,47] that the broad 
component of the signal for water diminishes with decrease in the distance 
between the silica particles. It has been suggested that it is due to the formation 
of amorphous modification of ice near the surface. 
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Fig. 4. Dependencies of 'H NMR spectra of water in frozen aqueous suspensions of 
fumed silica A, on (a) temperature and (b) silica concentration: 1 % w/w (1); 2 96 w/w 
(2); 6 96 w/w (3); and 20 96 w/w (4) 


Figures 5 and 6 show the AG(d) and F(x) graphs for silicas differently pre- 
treated. These graphs reflect the type of the radial functions relating to changes 
in the free energy of water in the adsorption layer. Two segments may be 
revealed on the dependencies obtained: the segment of a quick decrease in the 
thickness of the unfrozen water layer in a narrow range of the AG value 
changes (temperature is about 273 K) and the one, where d decreases relatively 
weakly in a broad range of changes in the AG values. As was shown in section 
2, there are weakly and strongly bonded waters respectively. The total value of 
the free energy of the adsorbed water using linear dependences AG (d) and 
AG, (d) may be defined according to Eq. (5). The AG values calculated using 
data shown in Figure 6 are summarized in Table 1. 
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Fig. 5. Effect of the pretreatment procedure of fumed silica on the shape of radial 
function of variation in the free energy of adsorbed water 


180-4— 
160 


X (nm) 
Fig. 6. Effect of adsorbent pretreatments on the radial function of adhesion forces for 
aqueous suspensions of fumed silica 


Tab. 1. Influence of silica pretreatment on the characteristics of the interfacial water 


Note. Unfrozen water layer thickness d are shown in statistical monolayers. 
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The data given in Table 1 show that the close values of AG; and the 
thickness of the strongly bound water layer have been obtained for all the 
adsorbents with the exception of A,, For weakly bound water, the thickness of 
the layer of structurally ordered liquid is minimal for A, and maximal for A,. 
For A, and A, samples, the ds; AG, and y parameters have not been 
determined taking into account superposition of the narrow and wide 
components of the signals (Figure 5). The wide component of the bonded water 
signal disappears as a function of rise in the bulk density of silica. However, the 
maximal thickness of the water layer structurally ordered by the surface is 
observed on slight compaction of the material. Thus, for practically equal 
values of C, for A, and A,., the thickness of the bound water layer varies more 
than twice. The dependencies obtained may be explained by the fact that the 
characteristics of the water layers structurally ordered by the surface depend on 
interparticle interactions. Fumed silica particles in the aqueous suspensions are 
present in the form of aggregates (100-500 nm) and agglomerates (> 1 um) 
[48]. These swarms are not stable and even exposure to negligible external 
factors (e.g. shaking of suspensions) gives rise to dispersion as well to 
amalgamation of particles. The amorphous form of frozen water is likely to be 
formed only at the appropriate distances between the particles and favorable 
dimensions of the secondary particles of the adsorbent. Therefore, the influence 
of the adsorbent pretreatment on its adsorption characteristics should be taken 
into account while studying adsorption of biopolymers on dispersed materials. 
By means of modifying the pretreatment conditions for the adsorbent one may 
vary its adsorption capacity in a wide range of values. This will have an effect 
on the adsorption value of biological macromolecules. 

The shape of radial dependence of the adhesion forces shown in Figure 6 for 
fumed silica in the aqueous medium differs significantly from similar 
dependences measured directly as forces between two crossed cylinders [49- 
-51]. Clearly from these results (Figure 6), the adhesion force value decreases 
practically linearly with the distance (not far from the surface). At the same 
time, the interaction force between two cylinders changes in the sign at a short 
distance between them due to too much decrease in the liquid layer thickness 
and can be described as follows F = x^. When x is smaller than the adsorbed 
layer the work is made against the adsorption forces; and the structure of this 
layer is destroyed. In the experiments on freezing-out of the bulk phase, the 
adhesion force F = dG/dx can be determined as AG/x (where x is the distance 
between the surface and the outer boundary of the unfrozen water layer). While 
the energy of interaction between the oxide surface and adsorbed water 
molecules typically increases with decreasing number of molecules in the 
interfacial layer, change in the sign of F is not observed. One can assume that 
such a method (however, as an indirect one) to estimate the work of forces of 
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adhesion and F(x) is more adequate than that based on direct measurements 
between two cylinders, as a cylinder can change the structure of the interfacial 
layer near the surface of second one producing additional errors in the 
measurements. 

Adsorption of polymers, e.g., poly(vinyl pyrrolidone (PVP), on fumed silica 
changes the free energy of the interfacial water. The adsorption potential 
distribution computed in respect to the disturbed (unfrozen at T < 273K) 
interfacial water using the 'H NMR measurements with entire freezing-out of 
the bulk water and _ layer-freezing-out of the interfacial water at 
210K < T < 273K (linear portions of the AG(C,,,) graphs were extrapolated to 
the X axis at AG — 0 in Figure 7a to determine the amounts of strongly C’,,, 
and weakly C",, bound waters) shows that the concentration of weakly bound 
water (Figure 7b, A « 1 kJ/mol) are lower for the pure silica suspension or the 
pure PVP solution. The decline in a linear portion of the AG(C,,,) graph at small 
AG (Figure 7a) corresponding to weakly bound water is larger for pure silica 
suspension than that for PVP/silica; i.e., the boundary between disturbed and 
bulk waters is more robust without PVP. However, in the case of pure PVP 
solution, the amounts of weakly bound water cannot be estimated, while nearly 
linear AG(C,,,) allows one to calculate only the amount of the strongly bound 
water, as the PVP capability to disturb the water to a large distance is low. This 
effect can be also caused by formation of PVP oligomers in the pure aqueous 
solution, whose "outer surface area" is small. The f(A) intensity (Figure 7b) 
increases at Cpyp = 5 wt. 96, but at Cpyp = 1 wt. 96, the opposite effect is 
observed in respect to f(A) at A > 1 kJ/mol, C',, and the free surface energy Ys 
in comparison with those for the individual silica suspension (or C',, for 
PVP/physiological buffer solution (PBS)/silica). On the other hand, changes in 
AG as a function of the filled pore volume (channels in aggregates of primary 
particles and inter-aggregate space in agglomerates close to the outer surfaces 
of aggregates) (Figure 7c) and the pore radius (Figure 7d) show that the silica 
impact on water in pores at R, between 1 and 10 nm is significantly larger (as 
the free energy is lower) than that for the PVP/silica suspension (Cpyp = 1 wt. 
96, Y = Cpyp/Cs;o? = 0.17) due to reduction in long-range components of surface 
forces for silica particle swarms rearranged and shielded by PVP molecules. 
Clearly, in the case of large mesopores, unfrozen water can be located close to 
the silica particles, e.g., the outer surfaces of aggregates. These effects can be 
explained by formation of a denser PVP layer at the silica surfaces at low PVP 
concentrations (y « 0.2) resulting in the decrease in the amount of water 
strongly bound to the oxide surfaces. At large Cpyp = 5 wt. 96 (y = 0.83), 
a major portion of PVP is in the non-immobilized state or weakly interacts with 
the silica surfaces, but the amount of water weakly and strongly bound to PVP 
molecules increases significantly (Figure 7a). Thus, at y « 0.2, polymer 
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molecules adsorb strongly and the number of free pyrrolidone groups can be 
less than bonded ones (at y = 0.2, the relative number of free C=O groups in 
PVP adsorbed on silica from solution and then dried at room temperature is 
approximately 2/3, according to IR spectra); i.e., the adsorbed PVP layer is 
relatively dense, and oxide particles or PVP molecules can disturb only 
a relatively thin interfacial water layer. At great Cpyp, a significant portion of 
PVP molecules has free tails, which do not interact with the silica surfaces but 
effectively interact with water molecules, whose great portion is disturbed by 
polymer molecules. However, a decline in the AG dependence on C,,, V, or Rp 


(especially at AG » —0.7 kJ/mol) is small for a relatively thick layer of weakly 
bound water (Figure 7). 
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Fig. 7. (a, c, d) Changes in the free energy and (b) water adsorption potential 
distributions for pure PVP solution and suspensions of silica, silica/PVP, and 
silica/PVP/PBS computed on the basis of 'H NMR data; AG as a function of (c) filled 
pore volume and (d) pore radius for silica (6 wt. 96 A-200) and PVP/silica 
(1 wt. % PVP, 6 wt. % A-200) suspensions 
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4.2. Modified silica and mixed oxides. Mixed oxides are promised materials 
to utilize them as catalysts, polymer fillers, and adsorbents possessing surface 
Brónsted and Lewis acid sites [52-56]. Besides this, disperse particles of dust 
pollutants of industrial and natural origins are mostly the mixture of porous and 
non-porous oxides SiO», Al;O;, TiO», etc., containing carbon deposits on the 
surface, or chemisorbed molecules of organic compounds. In atmosphere, these 
particles serve as nuclei of condensation or crystallization for water vapour. 

In Figure 8, there are shown the dependencies of changes in the free energy 
versus the concentration of unfrozen water. The characteristics of the layers of 
water adsorbed on the surface of the studied samples and the free surface 
energy values calculated using Eq. (5) are represented in Table 2. 
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Fig. 8. Dependence of changes in the free energy of the interfacial water for fumed 
mixed oxides as a function of the unfrozen water concentration 
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It follows from Table 2 that for alumina-silica we have average value of the 
surface free energy between values corresponding individual oxides. For non- 
porous materials with a low surface concentration of alternate compounds 
(pyrocarbon on SiOz; SiO2/Al,03 on TiO;), we have a maximal value of the 
surface free energy. The most probable explanation of this fact is the 
polarization of the adsorbent surface caused by charge exchange between 
different phases. As was shown by Dukhin [12], upon the appearance of 
differently charged areas on the surface of disperse particles, there is 
a long-range component of the surface forces, which caused by the electrostatic 
field directed along the surface of particles. In this field, the dipoles of water 
molecules are under orientation effect. Remoteness of this type of interaction is 
determined by the values of charges and the distance between them. This 
causes en enhancement of the polarization of adsorbent particles and the 
appearance of the long-range component of the surface forces. In all the cases, 
only the surface polarization is the main factor, which affects the thickness of 
the adsorbed water layers. The increase in the carbon content to 40 wt.% 
reduces changes in the free surface energy of the adsorbent, however, it remains 
significantly higher than that observed for the parent silica. Additionally, the 
smaller the patches of the second phase densely grafted on the support, the 
larger the disorder in the interfacial water, which causes greater reduction in its 
free energy. 


Tab. 2. Characteristics of interfacial water on disperse oxides and pyrocarbon-silica 


Methyl silica lU 80 | | 
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| TiO; (le) * JP) 802 || — 12 uam — 341 / |i. 55.5612. | 
| SiO, +0.6%TiIO, | 95 | 
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Note. Unfrozen water layer thickness d are shown in statistical monolayers 
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Fig. 9. Radial function of adhesion forces for fumed mixed oxides 


It should be noted that in the case when the ys values are governed only by 
the hydrophilic properties of the surface, one can expect that the values of the 
free surface energy obtained for the carbosils should be lower than those 
obtained for the silica. In this connection, in the case of nonporous adsorbents, 
a dominant role in the formation of the structurized interfacial water layers 
relates to the tendency of the adsorbent surface to undergo polarization in an 
aqueous medium. For porous adsorbents, the distance between chargers on the 
surface is limited by the size of pores. So for this type of adsorbents 
polarization effects are insignificant and interaction with the aqueous medium 
is governed by the hydrophilic properties of the samples. The correlation 
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between the nanoparticle structure and the characteristics of water layers on the 
particle surfaces may help in estimation of the distances of dust pollution 
transfer and in predicting the effect of dust pollution on the precipitation level. 
Radial dependencies of adhesion forces for some modified and mixed oxides 
are presented in the Figure 9. 


5. MESOPOROUS ADSORBENTS 


5.1. Impact of hydrothermal treatment on the surface energy and adhesion 
forces. Silica gel Si-60 (Merck) (specific surface area Sper = 369 m/g, total 
pore volume V, = 0.753 cm’/g, average pore radius R, = 4.1 nm), Si-40 (Merck) 
(Sagr = 732 m'/g, V, = 0. 542 cm /g, R, = 1.48 nm) ded Si-100 (Merck) (Sger = 
332 m’/g, V =1. 153 cm’/g, R, = 6. 95 nm) were chosen as initial materials for 
hydrothermal treatment (HTT). at Tir = 200 °C resulting in samples Si-40-HTT 
€ 107 m Hg, V, = 0.497 cm /g, R, = 9.3 nm); Si-60-HTT (Sger = 47 m’/g, 

V, = 0.272 cm Mg, R, = 11.8 sm) and Si-100-HTT (Sger = 49 m’/g, 
V, = 0. 18 cm’/g, R, a2 vid) 

Structural features of hydrothermally treated silica gels Si-40, Si-60 and 
Si-100 were elucidated previously [57]. Vertical portions of the AG(C,,,) graphs 
(Figure 10a) are observed for the aqueous suspensions of initial silica gels 
Si-40 and Si-60 in contrast to non-porous adsorbents, initial Si-100 (Figure 10a) 
and all hydrothermally treated samples (Figure 10b) possessing larger pores. 
Appearance of this portion of AG(C,,) is due to unfrozen interfacial water in 
narrow pores with lowering temperature until T corresponding to a maximal 
AG value for this plot portion. On the other words, lowering T over some 
interval is not accompanied by water freezing, as this water is under strong 
action of the silica surfaces in narrow pores, when electrostatic fields of the 
opposite pore walls strongly overlap [58]. 

The concentration of strongly bound water C',, decreases with increasing 
pore size both for initial and hydrothermally treated silica gels (Table 3). The 
concentration of weakly bound water C",,, is larger for the HTT samples due to 
the enhancement of the pore size; however, the relationships between the C"uw 
values for different silica gels (initial and HTT) more complex than those for 
C,, (Table 3). Notice that the volume of unfrozen water (C",,* C',,) for 
HTT-200 samples is larger than their pore volume, maybe due to formation of 
thick water layers on the outer surfaces of treated samples. Also, filling of large 
pores by nitrogen could be incomplete at p/po = 0.98 (used to estimate V,) in 
contrast to filling by water in the aqueous suspensions. 
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Fig. 10. Changes in the Gibbs free energy of the interfacial water in the aqueous 
suspensions of (a) initial and (b) hydrothermally treated Si-40, Si-60, and Si-100; and 


corresponding (c, d) adhesion forces and (e, f) derivatives d(AG)/dC,,, as functions of 
the distance x (in À) from the silica surface 
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Tab. 3. Parameters of bound water layer for aqueous suspensions of silica gels, initial 
and after HTT at 200 °C 
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The free surface energy ys is maximal for Si-60, initial or HTT. This result 
can be caused by several reasons, as the free surface energy at the silica/water 
interfaces depends on both the nature (e.g. concentration of silanols, charge 
density) and the topography of the surfaces. HTT for all samples leads to 
increase in 3; but by different magnitudes and Ays (changes in ys due to HTT) is 
maximal (236 mJ/m?) for Si-100HTT (Table 3). Additionally, HTT results in 
a marked enhancement of AG,, which characterizes the free energy of the first 
monolayer of adsorbed water. This value for HTT samples is significantly 
greater than that for some other materials investigated previously, which can be 
due to a substantial increase in the concentration of silanols (=Si-O-Si= + HO 
— =SiOH + HOSi=) (concentration of silanols for different silicas in air can 
alter by several times [59]), capable to form strong hydrogen bonds with water 
molecules from the first layer at the larger number of such bonds per 
a molecule; i.e. the number of strongly bound water molecules in the first layer 
increases for HTT samples. However, for initial silica gels, the opposite result 
is observed. Notice that the F(x) function correlates with increasing pore size 
for HTT samples (Figure 10d), but such a correlation is absent for initial silica 
gels (Figure 10c). A similar effect is seen for the derivatives —d(AG)/dCyy 
versus x (Figures 10e and 10f). The -d(AG)/dC,, value is the analog of the 
chemical potential for isobaric freezing of water. Maxima of the dependence of 
-d(AGydC,, on x correspond to the distances from the silica surfaces 
characterizing more stable complexes (clusters) of the interfacial water. The 
position of these maxima shifts toward larger x with increasing pore size for 
both initial and HTT silicas (for HTT samples, the x values were computed in 
a model of a flat surface). These effects can be explained by the significant 
differences in the unfrozen water layers dependent on the pore size. It should be 
noted that for Si-40, the —d(AG)/dC,, maximum corresponds to a water 
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monolayer, but for Si-60 and Si-40HTT, this is 3-4 monolayers, and for 
Si-60-HTT and Si-100-HTT, it is five-six water monolayers (Figure 10). 


5.2. Influence of pyrocarbon on the adsorbed water layer structure. Silica 
gel Si-60 (Schuchardt Miinchen) (SG) was used to prepare carbon-mineral 
adsorbents by the pyrolysis of CH;Cl; in dynamic conditions using quartz 
rotary reactor at 823K for 0.5, 1, 2, 3, 4, and 6 h corresponding to different CS-i 
carbosils with various amounts of pyrocarbon Cc (Table 4) (synthesis of these 
samples was described in details elsewhere [60]). SG550 is Si-60 heated at 
550 *C for 24 h. 


Tab. 4. Structural and adsorption characteristics of adsorbents 


m?/g wt. %o n time, h cm /g 
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The characteristics of the water layers adsorbed on the carbosil and initial 
silica gel surfaces determined on the basis of dependences AG(C,,) are 
presented in Figure 11. Thickness of the bonded water layer (C,,") for the 
weakly bonded water and C, for the strongly bonded water) as well as free 
surface energy (Ys) compared in Table 5. 

It should be stressed that for porous materials, the amount of adsorbed 
(unfrozen) water cannot be much greater than the total volume of the adsorbent 
pores. Therefore, during a decrease in the total pore volume with increasing 
pyrolysis time, which results in formation of pyrocarbon deposits in the pores 
and on the outer surfaces of silica gel particles, there is a tendency to decrease 
the adsorbed water concentration while passing from SG samples to carbosils. 
If effectiveness of carbonization is small, then C,,," and AG,, can be determined 
quite adequately, but in case of other samples, they have a form of 
approximation. Appearance of inflections on the AG(C,,,) graphs can be caused 
by the heterogeneity of the adsorbent surfaces and limitation of the pore 
volume. 
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Fig. 11. Changes in the Gibbs free energy of the interfacial water in the aqueous 
suspensions of initial silica gel and mesoporous carbosils and radial function of adhesion 
forces for these adsorbents 


Hydration properties of the adsorbent surfaces at the adsorbent/water 
interface can be expressed the most completely through the surface free energy. 
It follows from the data in Table 5, that ys increases as a result of silica gel 
thermally treatment from 115 mJ/m^ to 150 mJ/m?. The initial stage of silica 
surface carbonization causes a decrease in the free surface energy whereby Y; 
decreases to 72 mJ/m? for CS-4. It is in agreement with the statement that the 
surface hydration properties are determined by the concentration of surface 
silanols, which can form the strong hydrogen bonds with water molecules. 
Adsorption sites for water on pyrocarbon deposits are oxidized carbon atoms, 
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whose concentration is much smaller than that of SiOH groups on the silica 
surface; therefore carbon deposits should weaken the hydrophilic properties of 
the hybrid adsorbents. However, it is stated that there is an increase in the 
adsorbent hydrophilic properties for the samples CS-5 and CS-6, which is 
evident in the increase of the values ¥ and C',,, and for CS-6, ¥ is even higher 
than that for the initial silica. 


Tab. 5. Characteristics of water adsorption layers from the 'H NMR data on the surface 
of mesoporous carbosils 


Adsorbent AG, C AG, 
kJ/mol mg/g kJ/mol 
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6. MICROPOROUS MATERIALS 


6.1. Activated carbons. Carbon adsorbent PS-0 produced utilizing plum stones 
(Wood Dry Distillation Works, Hajnówka, Poland) was used as the initial 
material. PS-O0 was subjected to hydrothermal modification (HTM) in 
a stainless steel autoclave (0.3 dm?). To prepare oxidized adsorbents PS-OX1 
and PS-OX2, 5 g of PS-0 was placed in a quartz thimble to the autoclave with 
20 cm? of 30 % H;O; and treated in the vapour phase at 523 K (PS-OX1) or 
623 K (PS-OX2) for 6 h. A reduced adsorbent (PS-H) was prepared on heating 
of PS-0 in a H; stream in a quartz flow reactor at 1073 K for 8 h [61]. 

Figure 12 shows variations in the free Gibbs energy as a function of the 
concentration of unfrozen water in the aqueous suspensions of the initial (curve 
1), reduced (curve 2), and oxidized carbons (curves 3 and 4). The AG(C,,) 
graphs exhibit the portions of the unfrozen water, which remain near constant 
over a wide AG range. This effect can be caused by water adsorption in 
micropores, as according to the theory of volume filling of pores [58], the free 
energy of water in micropores does not practically depend on its localization 
there. Adsorbed water can be frozen if its free energy becomes equal to that of 
ice at a given temperature. If such a temperature has not been reached, freezing 
of water does not occur, and the unfrozen water concentration remains 
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practically unchanged. The portion of the graph of AG as a function of C,,, to 
the left of the vertical line (Figure 12) corresponds to freezing of water in 
micropores; and that to the right is associated with water in larger pores or 
weakly bound water on the outer surface of the adsorbent particles. The points 
at the abscissa corresponding to the vertical lines determine the concentration 
of water adsorbed in micropores (Table 6, C,,,””). 


AG (kJ/mol) 


0 200 400 600 800 1000 
Unfrozen Water Conc. (mg/g) 


Fig. 12. Changes in the Gibbs free energy of the interfacial water in the aqueous 
suspensions of activated carbons 


Tab. 6. Characteristics of the interfacial water layers in the frozen aqueous suspensions 
of carbons 


Adsorbent AG max Louie: 
kJ/mol mg/g 


650 112 1201 


As is seen from the data in Table 6, the free surface energy (ys) and content 
of unfrozen water (C,,"^) decrease for PS-H. It is in agreement with the notion 
that such a treatment results in a decrease in the surface concentration of 
oxidized groups, which can form hydrogen-bonded clusters of water molecules. 
Additionally, for PS-H, a sharp decrease in the concentration of water adsorbed 
in micropores is seen. This effect can be explained by lowering of the 
accessibility of micropores for water molecules due to increase in the 
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hydrophobic properties of the adsorbent surfaces reduced by H;. When the 
carbon surfaces are treated with hydrogen peroxide (PS-OX1, PS-OX2), an 
increase in the bound water concentration is observed in comparison with PS-0 
due to formation of additional polar groups, which can form hydrogen-bonded 
complexes with water molecules in micropores and on the outer surface of 
adsorbent particles. 


6.2. Water adsorbed on zeolites. Zeolite samples [62,63] were synthesized 
according to the technique described elsewhere [64] on the basis of fumed 
alumina-silicas. Such a procedure allowed us to simplify the synthesis of 
zeolites while alumina-silica served as a source for formation of the both oxides 
in the zeolite crystal lattice. 

The zeolites were synthesized using alumina-silica with molar ratios 
SiO/ALO; equal to 131, 100, 26 and 19.5 for preparation of samples 1-4 
(Table 7), respectively. The template was tetrapropylammonium bromide. In 
order to prepare a reaction mixture, sodium hydroxide and template were 
dissolved in distilled water; then a suitable amount of alumina-silica was added 
to this solution mixed. In terms of oxides the component ratios in the reaction 
mixture were as follows: Na,/SiO, = 0.05, [(C4H;),N];O/SiO; = 0.04, 
H;O0/Na;0 = 300. 


Tab. 7. Structural characteristics of zeolites 


A X-ray diffraction study was made with a DRON-UMI (LOMO) 
diffractometer using Ni-filtered Cu Ka radiation. The comparison of the XRD 
patterns of the prepared samples with the literature data on ZSM-5 [65] shows 
that these samples correspond to ZSM-5, whose degrees of crystallinity were 
94-96 %. The specific surface area $ was measured by the method of thermal 
desorption of argon; the pore volume V, was determined using methanol 
adsorption isotherms. The values obtained are listed in Table 7. 

As follows from the data listed in Table 8, the concentration of water 
adsorbed in zeolite pores and the free surface energy of adsorbents (with the 
exception of the sample containing 3.15 wt. % of Al,03) decrease with the 
aluminum content. At the same time, the maximum decrease in the free energy 
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in a layer of bound water is practically the same for all the studied zeolites and 
about twice as large as the corresponding value for silicalite. Under these 
conditions, the total adsorbent pore volume is filled by water only for a sample 
at Canos = 1.96 wt.%. The observed regularity disagrees with the data presented 
in [66] for adsorption of water vapors in zeolites differing in the concentration 
of Al in their crystal lattice. The cause of the observed features may be two 
main factors, namely change in the structure of water polyassociates with 
increasing concentration of Al and variation of conditions for the competing 
sorption of water and gases in air. Really, in the situation of a high degree of 
hydration of the surfaces, water not only interacts directly with primary 
adsorption sites but also forms large clusters near Brónsted acid sites [67]. 
However, the high acidity is characteristic only of a relatively small portion of 
surface hydroxyl groups, which is evidenced for by the data of calorimetric 
titration [68]. With increasing concentration of aluminum, the content of such 
active OH groups decreases, and, as a consequence, one may also observe 
a decrease in the sizes of clusters of water molecules adsorbed on the surface. 
Thus, despite the increase in the concentration of primary sites for water 
adsorption, the enhancement of the concentration of aluminum can be 
accompanied with a reduction of the degree of hydration of the zeolite surfaces. 


Tab. 8. Characteristics of interfacial water layers in aqueous suspensions of zeolites 


If adsorption of water is affected in air, there may appear a second factor 
that affects the dependence of the degree of hydration of the surface on the 
aluminum content in a sample. Owing to the presence of air, the zeolite surface 
can adsorb simultaneously not only water but also molecules of nitrogen and 
oxygen that are not able to form hydrogen-bonded complexes with surface 
hydroxyl groups and are adsorbed predominantly on more hydrophobic portions 
of the surface. With increasing sizes of clusters of water adsorbed on the 
surface, some portion of the adsorbed gases must be removed from the surface. 
Then the value of adsorption of water should depend on the relationship 
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between the free energy of formation of water clusters and that of desorption of 
gas molecules from the surface. In this situation, the formation of water clusters 
on hydrophilic sites of the surface may be energetically less favourable than 
adsorption of gas molecules on hydrophobic portions of this surface. As the 
content of aluminum in the zeolite increases, the number of primary sites for 
adsorption of water increases too, but the water clusters formed nearby these 
sites make contact predominantly with hydrophobic portions of the surface 
created by siloxane bridges. As a result, one can observe an effect that 
manifests itself in a decrease in the free surface energy of zeolites with 
increasing content of aluminum in their lattice. 


7. CONCLUSIONS 


The free surface energy of unmodified and modified oxides in the aqueous 
medium may vary in the range of 100-1100 mJ/m? depending on the nature of 
the adsorbents, and the thickness of the interfacial water disturbed by the 
surfaces ranges from 3 to 45 statistical monolayers. This is connected with the 
fact that in aqueous medium a high contribution of polarization and structural 
components to measured values of surface energy is observed. The magnitude 
of these components depends not only on the morphology and chemical 
structure of adsorbent surfaces but also on the structure of inter-particle 
aggregates as well as on concentration of a solid phase in the aqueous 
suspension. For this reason for the same material, the thickness of the layer of 
adsorbed water may vary depending on conditions of preparation of the 
adsorbent in a wide range. The increase in the thickness of layer of surface- 
-structurized water leads to a decrease in the adsorption ability of the adsorbent. 

Maximum values of the free surface energy are observed for samples 
characterized by the heterogeneous surfaces (carbosil having a small amount of 
pyrocarbon deposits, the mixture of silica with methylsilica, amorphous silica 
previously wetted). Specially separated fragments of the surfaces of the above 
adsorbents are characterized by different values of dissociation constants of the 
surface hydroxyl groups and by different affinities to protons. At the same time 
a polarization of the adsorbent surface takes place. It can be assumed that the 
thickness of the layer of adsorbed water is strongly affected by the surface 
charge distribution nonuniformity connected with appearance of oppositely 
charged fragments, and the field created by these fragments influences the 
orientation of dipoles of water molecules. The appearance of surface charges is 
also possible in the zones of interparticle contacts. For this reason in the case of 
low dispersion materials, the values of the free surface energy and layer 
thickness of adsorbed water in the aqueous suspensions are not determined by 
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hydrophilic-hydrophobic properties of these materials but by capability of the 
surface to polarization. 

The maximum layer thickness for water adsorbed in porous adsorbents 
cannot substantially exceed their pore radius. Therefore, when the action radius 
for surface forces is smaller than the pore radius, the layer thickness for water 
adsorbed on the surface is governed mainly by the surface concentration of 
hydrophilic sites. If the surface potentials of the opposite walls of a pore are 
overlapped, all the water filling the pore is bound to the surfaces. The freezing 
temperature for water in such a pore is governed by the pore sizes and force of 
interaction between the surface and water. The free surface energy value 
measured for such adsorbents determines the intensity of interactions between 
water molecules and pore surface. In the case of chemically modified 
carbonaceous adsorbents the comparison between values of their free surface 
energies in an aqueous medium makes it possible to obtain valuable 
information on the constitution of their surface. 
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